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The pT-spectra for direct J /tp and ip' in hadroproduction at Tevatron energy have 
been calculated within the framework of the NRQCD formalism and the fragmenta- 
tion model in the collinear parton model as well as in the /cT-factorization approach. 
We have described the CDF data and obtained a good agreement between the pre- 
dictions of the parton model and the /cj'-factorization approach. We performed the 
calculations using the relevant leading order in Ug hard amplitudes and taking equal 
values of the long-distance matrix elements for the both models. 



I. INTRODUCTION 

During the years after the measuring of charmonium production! cross sections and po- 
larization effects for J / ip and ip' mesons at the Tevatron Collider |l| the phenomenology of a 
quarkonium production has a phase of intensive developments. Nowadays, it is understand 
that the heavy quarkonium production is a very complicated physical process needs many 
new theoretical ideas. Starting from the Color Singlet Model it was developed, so-called, 
the nonrelativistic QCD (NRQCD) formalism 0] to describe the nonperturbative transition 
of QQ-^iaii into a final heavy quarkonium. 

The perturbative fragmentation functions for partons which split into heavy quarkonia 
have been obtained 4] within the framework of the NRQCD formalism. It is supposed that 
the fragmentation model is more adequate for the description of a quarkonium production 
at the large quarkonium transverse momentum px ^ ^qq than the fusion model. 
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In the charmonium production at the energy range of the Tevatron Colhder we deal with 
the gluon distribution function from a proton which is taken at a very small x but the 
relevant virtuality /i^ ~ ^qq + Pr is large. In the region under consideration the collinear 
mrton^niodel can be generalized within the framework of the fcr— factorization approach 
This fact leads to some interesting effects in the quarkonium production at 
he hig h e nerg ies, which were discussed ten years ago in Ref. 0] and recently in Refs. 
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In this paper we calculate the py-spectra of the unpolarized direct J/ip and ip' mesons 
produced via the fragmentation mechanism. In case of direct J/ip and ip' it is supposed 
that the production via the gluon fragmentation into the color-octet QQ [^51,8] state is a 
dominant contribution Q]. We compare the predictions which are obtained as in the 
coUinear parton model as in the /cT^-factorization approach. In both cases we performed 
calculations with the hard amplitudes in the leading order of the QCD running constant a^. 
So, in the collinear parton model we take into consideration the partonic subprocess: 

9 + 9-^9 + 9- (1) 

In the /cy-factorization approach we take into consideration the subprocess with off-shell or 
reggeized initial gluons: 

9* + 9* -^9- (2) 

n 

We have described the CDF data 1] and obtained a good agreement between the parton 
model and the fc^-factorization calculations based on the hard amlitudes for subprocesses 
(1) and (2) with the equal values of the long-distance matrix elements < O'^^'^fSi, 8] > and 
< 0'^'[^Si,8] > in the fragmentation functions Dg^j/^^^'{z, ^'^) for the both models. The 
QCD evolution of the fragmentation function is described by the homogeneous equation 
with a boundary condition proportional to the delta function 6{1 — z). 

The results obtained in the paper within the fc^-factorization approach based on the 



fragmentation function model differ from results of Refs. 



11 



12, hj, which were obtained 



using the gluon fusion model. The possible reasons of a disagreement will be discuss below. 
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II. NRQCD FORMALISM 

Within the framework of the NRQCD, the cross section or the fragmentation function 
for a quarkonium H production can be expressed as a sum of terms, which are factorized 
into a short-distance coefficient and a long-distance matrix element 

da{H) = Y.da{QQ[n]) < >, (3) 

n 

D{a -^H) = J2D{a^ QQ[n]) < > . (4) 

n 

Here the n denotes the set of color and angular momentum numbers of the QQ pair, which 
production cross section is da{QQ[n]) or which fragmentation function is D{a QQ[n]). 
Last ones can be calculated perturbatively in the strong coupling Of course, in the case 
of production in a hadron collision, the short-distance cross-section da have to be convoluted 
with the parton distribution function from the hadrons. The nonperturbative transition from 
the QQ state n into the final quarkonium H is described by a long-distance matrix element 
< (9^[r7,] > which have to be calculated using non perturbative methods or determined 
from experimental data. The fit of the Tevatron data for the p^-spectra of J/ip, ip' and 
charmonium states have been done recently by different authors (see review in Ref. 



As it was shown in Refs. 



igLAc 



15|, at the large charmonium transverse momentum [p^ > 7 



GeV) the gluon fragmentation into color-octet state 

g*^QQ[^Su8] (5) 

gives dominated contribution to the direct J/ip and ip' hadroproduction. The values of the 
color-octet matrix elements under consideration are following < (9"^/'^[^S'i, 8] >=4.4 ■ 10~^ 
GeV^ < 0^'[^Si,8] >=4.2 ■ lO^^ GeV^ Q- Note, that the fit of CDF data for the pr- 
spectrum of a direct J/ ip production within the framework of the fusion model in the collinear 
parton model gives the another numerical value of the long-distance matrix element < 
OJ/^^^Si,8] >=1.2 ■ 10~^ GeV^ 0. It will be important to compare the results obtained 
in the /Ey-factorization approach based on the fusion and fragmentation models. 



III. FRAGMENTATION FUNCTION 



The gluon fragmentation into ^Si charmonium state is determined by the color-singlet 
(Fig.l,a) and the color-octet (Fig.l,b) contributions. The previous analysis has shown that 
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the probability of a gluon fragmentation into the color-singlet state is only a small part of 
the probability of a fragmentation into the color-octet state [143. '^^^ leading order in the 
Us fragmentation functions for the transition (4) are known and can be written at the scale 

n 

/i^ = /ig = 4m^ as follows : 



where 



48m3 



6{l-z) 



(6) 
(7) 

(8) 



(9) 



It is obviously that the probability of the gluon fragmentation into the longitudinally 
polarized charmonium is negligibly small and J/ ip or tp' mesons have transverse polarization. 

The fragmentation functions (7) and (8) are evolved in /z^ using the standard homoge- 
neous DGLAP evolution equation [17.] : 



J z X Z 
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(10) 



where Pgg{x) is the usual leading order gluon-gluon splitting function. To solve equation 
(9) we use the well known method based on Mellin transform. It is easy to obtain that the 
Mellin-momentum at the scale /i^ can be written as follows 



Dg{n, fi^) = Dg{n, fil) exp 



2tt 7^*2 ^: 



-asifx^) 



where 



(11) 



(12) 



6 n{n-l) {n + l){n + 2) 

n 1 

^i(n)=E-- 

j=iJ 

In the one-loop approximation for the running constant a;s(yU^) with the three active flavors 
(6o = 9) one has 



Air 



6olog(/iVA2) 



(13) 
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and equation (10) can be presented as follows 

Dg{n, /i^) = Dg{n, /Iq) exp 
Taking into consideration that 



bo ' "Vl0g(/i2/A2 



(14) 



^J(-'/^o) = ^^^<(^^/^'^'p5„8]> (15) 
we have performed inverse Mellin transform numerically using following rule 

Dg{z, /x^) ^ Bl{z, 11^) = 1^ dnz-'^Djin, /x^). (16) 



1 r 1 

Dg{z,f/) = - dt Im e*^;z-"Dj(n,/i2) , (17) 
n Jo L J 



The integration contour C can be transformed such as 

1 /""o 

TT Jo 

7T 

where n = c + te^'^ and c ~ 2, = — . In Fig. 2 the obtained fragmentation function 
Dg^j/^{z,^i^) multiplied by 10"^ is shown at the different /i^ = 30,100 and 300 GeV^. We 
see that our result at the //^ = 300 GeV^ agrees well with the same one obtained in Ref. Q. 
In a stage of convoluting of the fragmentation function Dg^j/^{z, fi^) with the partonic cross 
section for the subprocesses (1) or (2) we will use the following definition for the variable z: 

(18) 



2Eg 

Thus we consider that the massless parton fragments into the massive meson. As it will be 
seen the definition (18) is more correct at a not so large px than the massless one between 
the parton 4-momentum and the meson 4-momentum 

We suggest also that meson has a small transverse momentum respectively initial gluon jet 
and we approximately can accept that in the laboratory frame 

= Vg, (20) 

where ri^,r]g are the meson and gluon pseudorapidities. 

Within the framework of the fragmentation model, the meson production cross section 
and the relevant gluon production cross section are connected as follows 

S'igg J/^X) = J dzDg^j/^{z,jj,'^)a{gg gg) (21) 



or 



a{g*g* ^ J/^X) = j dzDg_^j/4z, fi^)a{g*g* ^ g). (22) 
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IV. LEADING ORDER HARD AMPLITUDES 

The squared amplitude for the partonic process (1) is well known and it can be presented 
as follows 

\M{gg ^ gg)\^ = 1871 ^j^^ , (23) 

where s, i and u are usual Mendelstam variables. 

There two approaches for calciilation of a partonic amplitude for the subprocess (2) in 

The effective Feynman rules for processes with off- 



the /cT-factorization approach 
shell gluons were suggested in Ref. 
polarization 4-vector as follows 



The special trick is a choice of the initial gluon 



In Ref. 



e^ikr) = (24) 

\Kt\ 

the initial gluons are considered as reggeons (reggeized gluons) and the effective 
Reggeon-Reggeon- Gluon vertex function was obtained 

^-iM- . ( Jry + 2|^) - , 2^) , (25) 

where Pi = ^(1,0,0,1) and Pi = ^(1,0,0,-1) are the colliding protons 4-momenta, 
ki = xiPi + kiT and k2 = X2P2 + k2T are the initial gluons 4-momenta, kx = (0, A;t,0), 
k = ki + k2 is the final real gluon 4-momentum. It is easy to show that vertex function 
C^{ki, k2) satisfies the gauge invariance condition {ki + k2)\C^{ki, k2) = 0. 

Omitting the color factor f'"^ we can write the amplitude of the subprocess (2) accordingly 
Ref. as follows 

M = -ge^{k)-^^^^- [{k + A;i),(7am + (-^^i + k2)xg^u + {-ki - k)^g,x] ■ (26) 
We have obtained after simple transformations that 

M = f^^^^ xiX2sCxiku /C2), (27) 

2\kiT\\k2T\ 

where 

opA 2P^ 

C\ki, k2) = -{ki - k2f + ^(kl + X1X2S) ^ (kl + X1X2S) = C\ki, k2). (28) 



vertex function and amplitude: 



C\ki,k2)Cx{kiM) = (29) 

XIX2S 
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and 



\M{g*g*^g)\'^ = -'Ka,fT. (30) 



where = {kix + A;2t)^ = X1X2S, px is the transverse momentum of the final gluon. 

V. CROSS SECTIONS FOR THE PROCESS pp J/ij{7p')X 

In the conventional collinear parton model it is suggested that hadronic cross section, in 
our case, a{pp J/iIjX,s), and the relevant partonic cross section ^{gg J/tpX^s) are 
connected as follows 



{pp ^ J/ijj^s) = J dxi J dx2G{xi, i/)G{x2, n'^)a{gg J/^Ij,s), (31) 



where s = X1X2S, G{x, p^) is the collinear gluon distribution function in a proton, xi^2 are the 
fraction of a proton momentum, p? is the typical scale of a hard process. The evolution 
of the gluon distribution G{x,pL^) is described by DGLAP evolution equation [l7|. 

In the /c^^-factorization approach hadronic and partonic cross sections are related by the 
following condition 0, Q : 

a^^(pp J/^X) = / ^ / dklr J ^^xu kf^, p') (32) 



where d'{g*g* — > J/ipX,s) is the J/tp production cross section on off-shell gluons, kl = 
k\rp = —k\rp^ k\ = k^rp = — k^rp , s = X ix 2 s " (/cij" + /c2t)^; 01,2 are the azimuthal angles in 
the transverse XOY plane between vectors kiT{k2T) and the fixed OX axis {p^ G XOZ). 
The unintegrated gluon distribution function $(xi, A;^^, /i^) satisfies the BFKL evolution 



equation 



Our calculation in the parton model is down using the GRV LO f2Q[ and CTEQ5L 
21[ parameterizations for a collinear gluon distribution function G{x^p^). In case of the 



/c^^-factorization approach we use the following parameterizations for an unintegrated gluon 



distribution function ^(xi, kij,^ u _)- JB by Bluemlein |22|, JS by Jung and Salam |23|,KMR 
by Kimber, Martin and Ryskin [2^. The direct comparison between different parameteri- 
zations as functions of x, k^ and /i^ was presented in paper 
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The doubly differential cross sections for the process pp J/ip{ip')X can written as 
follows 



— — = dxi dzG[xi,n )G[X2,H )Dg_,^{z,fi ) " ^ j — ■ r, (33) 

where 

U = -Vs{Eg - Pg^), t = -^/s{Eg + Pg^), X^ = 7^^—, Xl^miu = 

U + XiS s + t 

i = xit, u = X2U, s = X1X2S; 
d^""^ f.D^.. t2 „2^^,. t2 „2^n „2^ Eg \M{g*g*^g)\^ 



j dz j d(f)i J dP^T^{xi,klT^ i/)^{x2,klrp, i/)Dg^^{z, jji'^ 



dr]^dp^T J J J PstE^ X1X2S 

(34) 



where 



r - r Eg+ Pgz Eg - 

k2T = PgT - kiT, Xi = , X2 " 



The energy (Eg) of a fragmenting gluon and the energy of a final meson [E^) are related as 
follows 

eo,«,,, ..^2a..,(e.p,-,,),. 

VI. THE RESULTS 

We compare our predictions with the CDF data [Ij for unpolarized direct J/ip and ip' 
production. The direct J/ip cross section does not include contributions from B meson 
decays into J/ip as well as from Xc radiative decays. For direct ip' production, indirect 
contribution from B meson decays are removed. Our results obtained in the coUinear parton 
model do not depend on a choice of the parameterization for the gluon distribution function 



and coincide approximately (±10-20%) to the results obtained in Ref. jl5| using a similar 
approch. We can see in Figs. 3 and 4 that the curves denoted as "GRV", which were 
obtained in the collinear parton model, are below the data especially at the small pt, were 
the gluon fusion into ^Sq and ^Pj color octet states is dominant [l5. Ilfi|. 

The curves which were obtained in the fc-r-factorization approach strongly depend on a 
choice of the unintegrated gluon distribution function. In the region of a large pt, where 
the fragmentation approach is more adequate, the results obtained with JB j22] and KMR 
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2J| parameterizations coincide well. However, JS parameterization 23| predicts the val- 
ues smaller by factor 2, which are near the values obtained in the coUinear parton model 
with the GRV gluon distribution function. Such a way, the px spectra of direct J/ip and 
ip' mesons obtained in the collinear parton model and in the fcr-factorization approach are 
approximately coincide. Note, that our conclusions disagree with the previous results ob- 



tained in Refs. [11 
data accordingly 



U, 12, 3 



_using the fusion model. Opposite our result, the fit of the CDF 
needs strong suppression (in 10-30 times) for the long-distance 
matrix elements < (^■^/'^''^'['^S'l, 8] > to compare the values obtained in the collinear parton 
model. 

There are several reasons of a such disagreement. At first, we used the fragmentation 



functions, which take into account effectively 
lution equation. At second, in Refs. |ll|, lia . 



ligh order corrections via the DGLAP evo- 
13| the argument /x^ of the strong coupling 



constant Q;s(/i^) is equal to klj, or fc n^, our choice is /i 



JJrp 



+ 4m^. This fact gives the 



additional factor 3 in a cross section [13[. At third, in Refs.|lll. Il2| the KMS|25| parameter- 
ization for an unintegrated gluon distribution function was used. We have shown (Figs. 2 
and 3) that the difference between predictions based on the different parameterizations may 
be about factor 2-3. 

The obtained results for the direct unpolarized Jjip and ip' hadropro duct ion as well as 
our previous results for the J/ip photoproduction at HERA energies show that the 
predictions for spectra on with the LO in «<, hard amplitudes in the /cT^-factorization 
approach coincide well to the predictions with the NLO in as hard amplitudes in the collinear 
parton model. The NLO in as calculation for the J/ip photoproduction cross section was 
performed in Ref. 

It is obviously that the NLO hard subprocess for a gluon production in the 
/ct— factorization approach with off-shell initial gluons is the LO subprocess used in the 
collinear parton model with on-shell initial gluons: 



g +g^g + g■ 



An amplitude of the subprocess (34) has infrared singularities even at the large value of pt 
for the final gluon which fragments into a meson. Opposite, in the collinear parton model 
the both gluons to be hard in similar case. The procedure of a calculation of an aniplitude 
for the subprocess (34) in the /cy-factorization approach has been suggested in Ref. jsf , were 
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initial gluons are considered as reggeons and the infrared divergencies are removed. It should 
be interesting to calculate J/ipltp') production cross section using the results of Ref. j^. This 
study is in progress. 
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FIG. 1: Diagrams used for description LO in the ag fragmentation g* J/ipgg (color-singlet 
state, A) and g* ^ J/ip (color-octet state, B). 




FIG. 2: Dg^j/^{z,iJ.^) at the /x2=30, 100 and 300 GeV^ as a function of z. 
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FIG. 3: The spectrum of direct J/V' on pT at the \/5=1800 GeV and \r]\ < 0.6. The data points 
are from [1]. The B is the J/V' lepton branching. 




FIG. 4: The spectrum of direct ij)' on pT at the \/S'=1800 GeV and |r/| < 0.6. The data points 
are from [1]. The B is the if)' lepton branching. 



